Spin transport and quasi 2D architectures for donor-based quantum computing 
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Through the introduction of a new electron spin transport mechanism, a 2D donor electron spin 
quantum computer architecture is proposed. This design addresses major technical issues in the 
original Kane design, including spatial oscillations in the exchange coupling strength and cross-talk 
in gate control. It is also expected that the introduction of a degree of non-locality in qubit gates 
will significantly improve the scaling fault-tolerant threshold over the nearest-neighbour linear array. 

PACS numbers: 03.67.Lx 
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The Kane paradigm of donor nuclear spin quantum 
computing in silicon (J, based on single atom placement 



coherent electron spin transport 



fabrication techniques [j, LJ|, is an important realization 
of Feynman's original concept of nanotechnology in the 
solid-state. Variations on this theme include electron 
spin qubits (j. and charge qubits [7]. There are sig- 
nificant advantages of the donor spin as a qubit, including 
uniformity of the confinement potential and high number 
of gate operations possible within the electron ^pin co- 
herence time, measured to be in excess of 60ms yj. Con- 
sequently, there is great interest in donor-based architec- 
tures and progress towards their fabrication |fi ll(i lll|. 

It is often assumed that solid-state designs should 
be inherently scalable given the capabilities of semi- 
conductor device fabrication. In reality this weak- 
scalability argument should be replaced with a stronger 
version as scalability of a given architecture is consid- 
erably more complex than fabricating many interacting 
qubits. Fault-tolerant scale-up requires quantum error 
correction over concatenated logical qubits with all the 
attendant ancillas, syndrome measurements, and classi- 
cal feed-forward processing. Both parallelism and com- 
munication must be optimised [12]. Only by considering 
such systems-level issues in conjunction with the underly- 
ing qubit physics will the requirements of quantum com- 
putation in a given implementation be understood, and 
new concepts generated. In this paper we introduce a 
new mechanism for coherent donor electron spin state 
transport, and in a similar design path to the QCCD ion 
trap proposal [13], we construct a 2D donor architecture 
based on distinct qubit storage and interaction regions. 

The significant interest in scaling up the donor-based 
solid-state designs, has led to a number of works consid- 
ering these scalability issues. As a result, several serious 
problems have been identified, including: sensitivity of 
the exchange interaction and control to qubit placement 
(at the 2-3 lattice site level) 0,0, [HI, qubit control and 
fabrication limitations associated with high gate densi- 
ties ll7L fl8H . spin readout based on spin-charge transduc- 
tion |l|,ll2|, and the communication bottlenecks for linear 
nearest neighbour (LNN) qubit arrays 

The issue of local versus non-local fault-tolerant oper- 
ation is non-trivial I2J A recent surprising result 
is that Shor's algorithm can be implemented on a LNN 
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Figure 1: Top view of the 2D donor electron spin quantum 
architecture for the case of Si:P, incorporating coherent trans- 
port by adiabatic passage (CTAP). 



circuit for the minimal qubit case with no increase at 
leading order in the circuit gate count or depth (23I E3|. 
However, at the systems level one expects a linear near- 
est neighbour qubit array to suffer from swap gate over- 
heads, particularly when concatenated qubit encoding is 
employed. The general analysis in [22] shows that local- 
ity forces the threshold down inversely with the physical 
encoding scale. Recently, the extent of the LNN penalty 
has been estimated to bring the threshold down by two 
orders of magnitude compared to the non-local case |25l |. 

For the Kane, or related donor based architectures, all 
of the above implies the imperative of finding ways of 
traversing the linear array constraints, as the most effec- 
tive way to improve the threshold and tackle the technical 
problems listed. An important step in this direction is 
the proposal for sub-interfacial transport of electrons in 
a one dimensional array jifj. This design has many de- 
sirable features, digitising the single and two qubit gate 
problems in an elegant way, but also has problems with 
scalability due to the relative closeness of gates (l^ |. 
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Figure 2: Top: Schematic of the one-electron triple donor 
system 3D 2+ based on P donors in silicon. Two of the donors 
are assumed ionized, the other neutral. Bottom: multi-donor 
CTAPn straddling schemes. 



The 2D architecture introduced here requires relatively 
low gate densities and specifically address the problems 
listed above. In Fig. [J] the geometry is shown for the spe- 
cific case of the exchange-interaction based Kane archi- 
tecture. We note that the transport ideas presented here 
allow for a similar, but non-trivial development for the 
digital-Kane case. A buried array of ionised donors pro- 
vide pathways for coherent transport of electron spins for 
in-plane horizontal and vertical shuttling (dashed-border 
sections) of qubit states into and out of the interaction 
zone. The overall gate density is low compared to the 
Kane case, and can be further reduced by increasing the 
transport pathway length (Fig. 0). Initially all gates in- 
hibit tunnelling along any given channel. Coherent spin 
transport along one segment is achieved by adiabatically 
lowering the barriers in a well defined sequence to effect 
coherent transfer by adiabatic passage (CTAP) without 
populating the intervening channel donors |2jj. We show 
that with appropriate donor separations, the shuttling 
time can be in the nanosecond range for one section. In 
Fig. [I] the coherent transport scheme is defined for the 
minimum number of donors. Higher order schemes with 
more donors reduces the gate density (see Fig. |2j). 

Logic gates are carried out in interaction zones distinct 
from qubit storage regions - shown in Fig. G]are the can- 
nonical A and J gates for electron spin based qubit con- 
trol at the microsecond level 0|. After mandatory preci- 
sion characterisation j2^| [2<j, interaction regions with un- 
acceptably low couplings can be identified and bypassed in 
the circuit flow, thereby avoiding bottleneck issues arising 
from the sensitivity of the exchange interaction to donor 
placement. This design allows for new variations on the 
theme, e.g. digitisation of hyperfine control [26], or intro- 
duction of local buried B-field antennae structures j^j, 
and space for SET readout techniques 0, [HI 0] . 



A schematic of the minimal three donor transport 
pathway is given in Fig. The triple- well system |1<j), 
|2<j), \3a) (a =|,|) facilitates coherent state transport 
from a\l j) + |) to a\3 |) + /?|3 |) without populat- 
ing the \2a) states. Techniques for coherent transfer by 
adiabatic passage are well known j^, an d for the donor 
system was proposed in (2g for the case of charge trans- 
fer. A superconductin g ve rsion of the three state case 
has also been proposed (23 . The system is controlled by 
shift gates, S, which can modify the energy levels of the 
end donors, and barrier gates, B^i+i which control the 
tunnelling rate between donors i and i + 1. 

Although the scheme we introduce here necessarily in- 
cludes spin, we first consider the zero field case and ignore 
spin degrees of freedom to illustrate the principles 
of CTAP in the one-electron three-donor system, 3D 2+ . 
The effective Hamiltonian for the 3D 2+ system is: 

H = A|2)(2| - ft (fii 2 |l>(2| + n 23 |2)(3| + h.c.) , (1) 

where Qij = flij(t) is the coherent tunnelling rate be- 
tween donors \i) and \j) and A = E^ — E\ = E<i — E%. 
The eigenstates of H (with energies £± and £q) are 
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sin 6i sin 2 11) + cos 62 12) + cos 61 sin 02 13), 
sin©i cos 62 1 1) — sin 62 1 2) + cosOi cos 62 1 3), 
cos8i|l) -sinOi|3), (2) 



where we have introduced 9i = arctan (^12/^23) and 
B 2 = arctan[2ftV / (^i 2 ) 2 + (^ 2 s) 2 /A]/2. Transfer from 
state |1) to 1 3) is achieved by maintaining the system in 
state I Do) an d changing the characteristics of \T>q) adia- 
batically (\£q-£±\ > |(i> |£>±)|) from |1) at t = to |3) 
at t = t max by appropriate control of the tunnelling rates, 
without population leakage into the other eigenstates. 

For the case of coherent spin transport we write the 
3D 2+ Hamiltonian in terms of spin/site operators as: 
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« = EE £-4^+ E M*hW (3) 



i=l cr = T, 



<ij> cr = T, 



and numerically solve for the density matrix, p(t), 
in the presence of a (dominant) charge dephas- 
ing rate T, assumed to act equally on all coher- 
ences. Without attempting to fully optimize con- 
trol we apply Gaussian pulses of the form ttij(t) = 
^•^ ax exp [—(t — Uj) 2 / (2wij)] , where Uj and Wij are the 
peak time and width of the control pulse modulating the 
tunnelling rate between position states \i) and \j). To 
simplify matters for initial simulations we set the max- 
imum tunnelling rates and standard deviations for each 
transition to be equal, i.e. £7™ ax = £7 max and = w, 
and set A = (these conditions can be relaxed with 
no effect on the conclusions of this paper). Transfer is 
then optimized when the width of the pulses equals the 
time delay between the pulses [34]. With total pulse time 
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Figure 3: Numerical simulation of the CTAP pulse scheme 
applied to a spin superposition at donor 1 at t = 0, demon- 
strating coherent transfer to the 3rd donor at t = t m ax- 

W 5 we choose w = tmax/ 8 so that t 12 = (t max + w)/2 
and ^23 = (tmax — w)/2. This ordering, where ^23 is ap- 
plied before Q12 is known as the counter-intuitive pulse 
sequence and has significant advantages in improving 
transfer fidelity over other pulse sequences |2g. In Fig.0 
we present results showing transport using the counter- 
intuitive pulse ordering for a spin superposition (phases 
relative to the untransported state). 

Generally, when the adiabaticity criterion is satisfied 
and the transport time is at least an order of mag- 
nitude faster than charge dephasing, the transport fi- 
delity is high. These results are consistent with those 
of Ivanov et al j^EI who considered the role of dephas- 
ing in three-state Stimulated Raman Adiabatic Passage 
(STIRAP). Although these competing timescales are es- 
sentially unmeasured at present, estimates |36l l37j for 
the P-P + charge dephasing time are of order 10ns and a 
value of 220ns was reported recently for a Si:P double-dot 
whereas sub-nanosecond tunnelling times are possi- 
ble due to the strong confining potential of donor nuclei. 
The CTAP transport time will be defined primarily by 
the gate-assisted tunnelling rate, which we calculate as 
follows. Using the TCAD package we compute the po- 
tential due to a surface B-gate bias and determine the 
donor electron wave function in an effective mass ba- 
sis, e.g. F±;' m (r) 

— ^n^mf^)!/^^)) about the six band 
minima where the ipn,i,m are hydrogenic orbitals with 
Bohr radius a± 7 and 7 = a±/a\\. Diagonalising the total 
Hamiltonian of the system, using pseudopotentials to de- 
scribe the silicon bandstructure, we obtain a generalised 
Kohn-Luttinger wave function: 
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<P(r,V)= ^c„, i)m (y)^F;^(r)e' k - r ^(r), (4) 

n,l,m [i=l 

where the Bloch states are u\^^(r) = A^^ (G)e* G " k/x . 
We form bonding and anti-bonding states \P±(r, V) = 
j/V(^L(r, V) ± ^(r, V)), normalised by A/", and compute 
the gap as shown in Fig. 0| for basis sizes 55 and 140 
(n max = 5 and 7). Comparison of the non-linear regions 
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Figure 4: Main: Energy gap for the P-P + system as a function 
of B-gate bias V& for R=30nm (depth 30nm below interface, 
lOnm gate width, basis sizes N = 55 and 140). Lower right: 
response of the P-P + inter-donor potential profile to the bar- 
rier gate bias V& = (0, ±500) mV. 



indicates that the range of validity is |V5>| < 200 mV. 

In contrast to what one expects for an isolated P- 
P + system in vacuum where the nodal structure of the 
bonding and anti-bonding states is simple, the non-trivial 
nodal properties of the donor electron wave function and 
the proximity of the oxide interface complicates the tun- 
nelling control. These calculations directly extend simi- 
lar effects noted in the ungated P-P + system j^j. From 
Fig. E| we see that for this configuration the tunnelling 
rate can be varied from zero at +100mV to ~ 10 GHz at 
-200mV, giving a gate assisted tunnelling time of 60 ps. 

Based on this value, CTAP simulations for 5, 7 and 
9 donor chains are presented in Fig. The adiabatic 
nature of the transport scheme provides an inherent ro- 
bustness, as evidenced in Fig. which shows a remark- 
able uniformity in the response to charge dephasing for 
the different path lengths once the adiabatic regime is 
reached. Another consequence is that inevitable varia- 
tions in tunnelling rates due to donor placement will 
not affect the viability of the scheme, as further simula- 
tions have explicitly verified. The extent to which T con- 
trols the transport fidelity is also clear, although we note 
that there is room for improvement through optimisation 
of control pulses and minimisation of charge fluctuations 
through fabrication development. Non-zero transport er- 
rors may require monitoring mechanisms for heralding 
successful transport, or an error correction protocol for 
transport loss. As intrinsic spin-orbit coupling for donor 
states in silicon is very low, dephasing of donor electron 
spin is dominated by spectral diffusion due to spin impu- 
rities and is mitigated by isotopic purification 01 . For 
the bound state spin-orbit coupling, at Vb ~ 200 mV 
we calculate from EqnQ the non-S components to be 
E n ,i>o,m \cn,i,m(V)\ 2 < 10~ 4 indicating that the devia- 
tion from the S sector is minimal. Together with the zero 
occupation of channel states, this suggests that charge 
dephasing will have a negligible second order effect on the 
spin coherence during transport. Decoupling of orbital 
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Figure 5: Transfer error as a function of charge dephasing 
rate and total transfer time for CTAP5 (solid line), CTAP7 
(long dashes) and CTAP9 (short dashes) for the case of 30nm 
end-donor spacings, and 20 nm between the central donors. 



and spin sectors has already given rise to demonstrations 
of coherent transport of electron spins over 100/im |4ll |. 

The basic layout of 2D donor arrays, with storage re- 
gions, vertical and horizontal transport pathways and in- 
teraction zones, allows us to explicitly consider designs 
for fault-tolerant operation. For example, we can arrange 
the logical qubit groups and ancillas so that the trans- 
port rails allow for non-local intralogical qubit interac- 
tions (qubit - ancillas) and LNN interlogical interactions. 
With inherent parallelism of operation, interlogical gates 
can then be applied transversally as required to imple- 
ment fault-tolerant gates. Another possibility with less 
stringent fabrication requirements is a linear qubit stor- 
age with transport and interaction rails either side. 

The optimum arrangement for fault-tolerant opera- 
tion requires sophisticated systems level simulations 01 
to determine the best use of this medium range quan- 
tum transport capability, and the corresponding improve- 
ments on the LNN threshold. In any case, it is clear 
that the introduction of coherent spin transport to donor 
quantum computing allows us to address many problems 
in the Kane concept, and consider scalable fault-tolerant 
architectures with low gate densities, room for SET struc- 
tures and control, and a bypass mechanism for low value 
exchange gates. One expects the realities of the silicon 
crystaline environment will necessitate the characterisa- 
tion of transport pathways, however, the precision re- 
quirements of the adiabatic CTAP mechanism would be 
far less than the quantum gate threshold. 
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